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1 DIS kinematics 



The event kinematics of deep inelastic scattering, DIS, are determined by the negative square of the four- 
momentum transfer at the lepton vertex, = —q^, and the Bjorken scahng variable, x = Q^/2P ■ q, where P 
is the four-momentum of the proton. In the quark parton model (QPM), the interacting quark from the proton 
carries the four-momentum xP. The variable y, the fractional energy transfer to the proton in its rest frame, 
is related to x and by y ~ Q'^/xs, where ^/s is the positron-proton centre of mass energy. Because the HI 
and ZEUS detectors are almost hermetic the kinematic variables x and can be reconstructed in a variety of 
ways using combinations of positron and hadronic system energies and angles iQ . 

Neutral current (NC) DIS occurs when an uncharged boson (7, is exchanged between the lepton and 
proton. In QPM there is a 1-f 1 parton configuration, fig. |l|a, which consists of a single struck quark and the pro- 
ton remnant, denoted by At HERA energies there are significant higher-order quantum chromodynamic 
(QCD) corrections: to leading order in the strong coupling constant, as, these are QCD-Compton scattering 
(QCDC), where a gluon is radiated by the scattered quark and boson-gluon- fusion (BGF), where the virtual 
boson and a gluon fuse to form a quark-antiquark pair. Both processes have 2-1-1 partons in the final state, as 
shown in fig. |l|. There also exists calculations for the higher, next-to-leading (NLO) processes. 

HERA provides a unique opportunity to study the scale behaviour of Quantum Chrmodynamics in a single 
experiment with a cleaner background environment than that obtainable at hadron colliders. The data presented 
here were taken in 1994 and onwards at the ep collider HERA using the HI and ZEUS detectors. During this 
period HERA operated with positrons of energy = 27.5 GeV and protons with energy 820 GeV. A detailed 
descriptions of the HI and ZEUS detectors can be foimd in refs. and respectively. 



2 Structure functions 

Perturbative QCD (pQCD) does not predict the absolute value of the parton densities within the proton but 
determines how they vary from a given input. For a given initial distribution at a particular scale, Qq, Altarelli- 
Parisi (DGLAP) evolution enables the distributions at higher to be determined. DGLAP evolution resums 
the leading log((5^) contributions associated with a chain of gluon emissions. At large enough positron-proton 
centre-of-mass energies there is a second large variable 1/x and, therefore, it is also necessary to resum the 
log(l/x) contributions. This is achieved by using the BFKL equation In the DGLAP parton evolution 
scheme the parton cascade follows a strong ordering in transverse momentum ^ Ptu-i ^ ■■■ ^ Pri^ 
while there is only a soft (kinematical) ordering for the fractional momentum Xn < Xn-i < ■■■ < xi (see 
figure ^.) By contrast, in the BFKL scheme the cascade follows a strong ordering in fractional momentum 
Xn <C Xn~i <C ... ^ 2^1, whilc there is no ordering in transverse momentum. 

At small X the dominant parton is the gluon and the description of the structure function is driven by the 
behaviour of the gluon. Because of gluon splitting, g qq, pQCD suggests the small x behaviour of the sea 
quark and gluon distributions are strongly correlated. 

The kinematic plane covered by HERA and the fixed target measurements is shown in fig. ^. HERA has 
increased the reach in by about 2 orders of magnitude and can also probe nearly 3 orders of magnitude 
further down in x. The low x region is correlated with low values of Q^. The differential NC DIS cross section 
is related to three structure functions: 

- lir(^+^2(:r, Q2) „ y2F^^^^ Q2) ^ Y^xF,{x, (1) 

where Fj- = 1 ± (1 — y)^. The structure function F2 in QPM is just the sum of the quark densities multiplied 
by the appropriate electric charge; F3 arises from the weak part of the cross section and is negligible for 

< 5000 GeV^, and Fj^ is the longitudinal structure function and only becomes important for y > 0.6. Hence 
by measuring the differential cross section at HERA one is effectively measuring the structure function F2 . 

The F2 measurements are shown in figs. ^ and |5| as a function of x and respectively. The error bars are 
at the 5-10% level and the normalisation uncertainty is ~2%. There is a steep rise of F2 with decreasing x in 
all bins, fig. ^. Scaling violations in are clearly seen in fig. ^. Both HI and ZEUS have performed next 
to leading order (NLO) QCD fits ||^, |^ based on the DGLAP evolution equations using both HERA and fixed 
target data. Fig. || shows that these QCD fits describe the F2 data well, though it should be noted that the 
data can also be satisfactorily described by the BFKL prediction 

The scaling violations from the HERA data allow an estimate of the gluon density xg{x) at low values of a;, 
whilst the fixed target data are used to constrain the high x region. The extracted gluon densities from the fits 
are shown in fig || for a fixed — 20 GeV^ . The error band shows the statistical and systematic uncertainty 
taking into account correlations and variations in the mass of the charm quark, 771^, and the strong coupling 



constant, a^. The results of the two HERA experiments are in good agreement and the extracted densities agree 
with the results of NMC for large x. The resulting gluon distributions show a clear rise with decreasing x 
and have a 15% uncertainty at a; ~ 5 x 10^^. These NLO QCD fits are also in good agreement with the global 
QCD analyses performed by MRS and CTEQ whilst the prediction from the dynamical evolution of 
GRV |l| is too steep for x < lO'^. 

3 Forward Jet Production in DIS 

The F2 measurements fail to distinguish between DGLAP and the BFKL approach to the QCD evolution. 
The hadronic final states are expected to give additional information. For events at low x, hadron production 
between the current jet and the proton remnant is expected to be sensitive to the effects of BFKL or DGLAP 
dynamics. A possible signature of BFKL dynamics is the behaviour of DIS events at low x which contain a jet 
that has a transverse momentum p\{j) ~ (so minimizing the phase space available for DGLAP evolution) 
and has longitudinal momentum fraction (of the proton) Xj^t that is large (in order to maximise the phase space 
for BFKL evolution) , fig ^. 

In Fig. 1^, recent data from HI jl^ and ZEUS [O are compared with BFKL predictions |jl^ and fixed order 
QCD predictions as calculated with the MEPJET program at NLO. The conditions ptU) — Q and Xjet ^ x 
are satisfied in the two experiments by slightly different selection cuts. HI selects events with a forward jet of 
PtU) > 3.5 GeV (in the angular region 7° < 9{j) < 20°) with 

0.5 < pt(j)VQ' < 2 , xjet ^ E,,t/Eproton > 0.035 ; (2) 

while ZEUS triggers on somewhat harder jets of Exij) > 5 GeV and r]{j) < 2.6 with 

0.5 < Pt{j)VQ^ < 2 , Xjet = Pz{j)/Eproton > 0.036 . (3) 

Fig. 1^ shows that both experiments observe a forward jet cross section which rises steeply with decreasing x 
with substantially more forward jet events than expected from NLO QCD (labelled as Born in fig. ^.) A BFKL 
calculation (the stars) gives a better agreement with the data. The overall normalisation in this calculation is 
uncertain and the agreement may be fortuitous. Indeed, it should also be noted that both experiments observe 
more centrally produced dijet events than predicted by the NLO QCD calculations. The ARIADNE Monte 
Carlo (CDM) model describes the steeply increasing jet cross section with decreasing x. The ARIADNE model 
does not have a strong ordering in transverse momentum in the QCD cascade, akin to BFKL type dynamics, 
although it does not make explicit use of the BFKL equation. Whilst those models that adhere to the DGLAP 
formalism (LEPTO and HERWIG) fail to predict this large growth. Further careful investigation is necessary 
before claiming that BFKL is the mechanism for this enhanced forward jet production. 



4 Fragmentation Functions 

Fragmentation functions represent the probability for a parton to fragment into a particular hadron carrying a 
certain fraction of the parton's energy and, like structure functions, cannot be calculated in perturbative QCD, 
but can be evolved from a starting distribution at a defined energy scale. If the fragmentation functions are 
combined with the cross sections for the inclusive production of each parton type in the given physical process, 
predictions can be made for the scaled momentum, Xp, spectra of final state hadrons. Small Xp fragmentation 
is significantly affected by the coherence (destructive interference) of soft gluons ijl^ , whilst scaling violation of 
the fragmentation function at large Xp allows a measurement of as [ p^ . 

A natural frame in which to study the dynamics of the hadronic final state in DIS is the Breit frame . In 
this frame the exchanged virtual boson is purely space-like with 3-momentum q = (0, 0, —Q), the incident quark 
carries momentum Q/2 in the positive Z direction, and the outgoing struck quark carries Q/2 in the negative 
Z direction. A final state particle has a 4-momentum p^ in this frame, and is assigned to the current region 
if p§ is negative, and to the target frame if p^ is positive. The advantage of this frame lies in the maximal 
separation of the outgoing parton from radiation associated with the incoming parton and the proton remnant, 
thus providing the optimal environment for the study of the fragmentation of the outgoing parton. 

In e+e~ annihilation the two quarks are produced with equal and opposite momenta, zt^/s/2. This can be 
compared with a quark struck from within the proton with outgoing momentum —Q/2 in the Breit frame. In the 
direction of the struck quark (the current fragmentation region) the particle momentum spectra, Xp = 2p^ /Q, 
are expected to have a dependence on Q similar to those observed in e+e~ annihilation |2^, ^ ^ at energy 



In fig ||the log(l/a;p) distributions for charged particles in the current fragmentation region of the Breit frame 
are shown as a function of . These distributions are approximately Gaussian in shape with mean charged 
multiplicity given by the integral of the distributions. As increases the multiplicity increases and the the 
peak of the distributions shifts to larger values of log(l/xp). Figure || shows this peak position, \o%(\j Xp)max, 
as a function of Q for the HERA data and of -^/s for the data. Over the range shown the peak moves from 

~ 1.5 to 3.3. The HERA data points are consistent with those from TASSO and TOPAZ and a clear agreement 
in the rate of growth of the HERA points with the e^e~ data at higher Q is observed. 

The increase of log(l/xp)maa; can be approximated phenomenologically by the straight line fit log(l/a;p)ma2; = 
61og(Q) + c also shown in figure 0. The values obtained from the fit to the ZEUS data are h = 0.69±0.01(stat)± 
0.03(sys) and c = 0.56 ±0.02^°;*°. The gradient extracted from the OPAL and TASSO data is 6 = 0.653 ±0.012 
(with c — 0.653 ± 0.047) which is consistent with the ZEUS result. This value is consistent with that published 
by OPAL, h ~ 0.637 ± 0.016, where the peak position was extracted using an alternative method ||2^. A 
consistent value of the gradient is therefore determined in DIS and e+e" annihilation experiments. 

Also shown is the statistical fit to the data when h ~ \ {c — 0.054 ± 0.012) which would be the case if the 
QCD cascade was of an incoherent nature, dominated by cylindrical phase space. The observed gradient is 
clearly inconsistent with 6=1 and therefore inconsistent with cylindrical phase space. 

The inclusive charged particle distribution, 1/atot da/dxp, in the current fragmentation region of the Breit 
frame are shown in bins of Xp and in fi g. HG . The increasingly steep fall-off, at fixed , towards higher 
values of Xp as increases, shown in figure |10|, corresponds to the production of more particles with a smaller 
fractional momentum, and is indicative of scaling violation in the fragmentation function. For > 80 GeV^ 
the distributions rise with at low Xp and fall-off at high Xp and high Q^. In figure |l^ the HERA data are 



compared at = s to e+e" data |25 , again divided by two to account for the production of both a q and q. 
In the range shown there is good agreement between the current region of the Breit frame in DIS and the 
e+e^ experiments. 

5 Summary 

This review gives a brief summary of a small sample of the QCD results coming from HERA. Charged particle 
distributions have been studied in the current region of the Breit frame over a wide range of Q^. These result 
show clear evidence for scaling violation in scaled momenta as a function of and supports the coherent 
nature of QCD cascades. The observed charged particle spectra are consistent with the universality of quark 
fragmentation in e~^e~ and DIS. 

The intriguing rise of F2 at small x can be well described using conventional DGLAP evolution equations. 
The data though can also be described by the BFKL approach thus giving rise to ambiguities how to treat 
QCD in this small x regime. In order to resolve these ambiguities jet production in the forward direction has 
been studied and the cross section for such jets is seen not to be reproduced by Monte Carlo models based on 
DGLAP parton shower evolution. 
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Figure 1: (a) QPM (b) QCDC and (c) BGF diagrams 




Figure 2: Parton ladder diagram in DIS. 
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Figure 3: The kinematical region covered by HERA and fixed target experiments 




Figure 4: Measurement of F2(x, Q^) as function of x. 
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Figure 5: Measurement of F2{x,Q^) as function of Q^. 
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Figure 6: The gluon density xg{x) at = 20 GeV^ extracted from NLO QCD fits by the Hf , ZEUS and NMC 
collaborations. 



a. 

^ 250 

f 

200 
150 
100 
50 

ft, 



H ' Preliminary 



7°<G,<20° 
E,>28-7GeV, Ptj>3-5GeV 
E,>llGeV, 160<G,<173° 
y>0.1, 0.5<PtVQ^<2.0 




ZEUS 1995 



05 0.1 0.15 0.2 0.25 0.3 0.35 

Cross Section, Sys and Stat errors added in quadrature X 1 



^ 120 
100 
80 
60 
40 
20 




10 



-3 



• ZEUS Data 
-|_ ARIADNE 4.08 
. LEPTO 6.5 (SCI) 
: HERWIG 5.9 



10 



-2 



(a) (b) 

Figure 7: Measurement of forward jet production as function of x by (a) HI and (b) ZEUS experiments. The 
errors due to the uncertainty in the jet energy scale are show in the shaded band. 
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Figure 8: Preliminary ZEUS results of the evolution of the 1/Ndn/ d\og{l/Xp) distributions with Q'^. Only 
statistical errors are shown. 



ZEUS 1994-97 PreUminary 

4 HQ) 





□ 


1 1 1 
OPAL (eV) 


1 1 1 1 1 1 






o 


TASSO (ee) 










TOPAZ (e*e") 










HI 1994 (eP) 








■ 


ZEUS 1994-97 Prelimiiiary 




* - 




























Fit to ln(1/xp)„„ = b InQ + c 










phase space (b=1) 





3 4 5 6 7 8 9 10 20 30 40 50 60 70 80 90100 

Q (GeV) 



Figure 9: Preliminary ZEUS results of evolution of the peak position \og{l/Xp)jnax with Q'^. 
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Figure 10: The inclusive charged particle distribution, 1/atot da/dxp, in the current fragmentation region of the Breit 
frame as measured by HI and ZEUS collaboratiouK. The open points represent data from e^e~ experiments divided by 
two to account for q and q production (also corrected for contributions from Kg and A.) 



